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This report covers the period from September 1, 1988 to January 31, 1989 for 
NASA grant NCC2-552 provided to ELORET Institute for support of Timothy J. 
Lee. During this period four projects were completed and each has resulted in a 
publication in a refereed journal. 

• The first project (in collaboration with G. E. Scuseria, A. C. Scheiner and H. F. 
Schaefer) involved the application of the high level single and double coupled cluster 
(CCSD) method to the atmospherically important ozone molecule. As discussed in 
the enclosed reprint, the CCSD method is the first single-reference-based electron 
correlation method to successfully reproduce the ordering of the uq and u> 3 harmonic 
vibrational frequencies. 

• The second project (in collaboration with C. W. Bauschlicher, S. R. Langhoff 
and P. R. Taylor) concerned the theoretical investigation of the barrier height for 
the chemical reaction F+H 2 -» FH+H. This reaction is one of the few for which 
accurate experimental information of the reaction dynamics exists and there has 
been some controversy concerning the barrier height. The enclosed reprint describes 
the highest level of theory yet applied to this system. 

• The third project was a completion of work started at the University of Cambridge 
in England during TJL’s postdoctoral studies. Much work has been done on neutral 
and cationic hydrogen-bonded complexes but very few high level calculations have 
been reported on anionic hydrogen-bonded complexes. The enclosed preprint (in 
press in the Journal of the American Chemical Society) describes a very detailed 
theoretical study of four anionic hydrogen-bonded complexes. 

• The fourth project (in collaboration with P. R. Taylor) represents a detailed study 
of the first generally reliable diagnostic for determining the quality of results that 
may be expected from single-reference-based electron correlation methods and will 
no doubt be very useful in the analysis of results in future theoretical studies. The 
enclosed preprint (in press in the International Journal of Quantum Chemistry, vol. 
S29) defines and reports our investigations of the diagnostic. 
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Abstract 

HF, H 2 0, CN“ and their hydrogen-bonded complexes have been studied using 
state-of-the-art ab initio quantum mechanical methods. A large Gaussian one- 
particle basis set consisting of triple zeta plus double polarization plus diffuse s and 
p functions (TZ2P - diffuse) was used. The theoretical methods employed include 
self-consistent-field, second-order Mpller-Plesset perturbation theory, singles and 
doubles configuration interaction theory and the singles and doubles coupled cluster 
approach. The FH-CN", FH ”NC” and H 2 0-CN", H 2 0-NC" pairs of complexes 
are found to be essentially isoenergetic. The first pair of complexes are predicted to 
be bound by ~ 24 kcal/mole and the latter pair bound by — 15 kcal/mole. The ab 
initio binding energies are in good agreement with the experimental values. The two 
pairs of complexes exhibit small structural differences with the N"H hydrogen bond 
being shorter than the analogous C "H hydrogen bond. The infrared (IR) spectra 
of the two pairs of complexes are also very similar, though a severe perturbation of 
the potential energy surface by proton exchange means that the accurate prediction 
of the band center of the most intense IR mode requires a high level of electronic 
structure theory as well as a complete treatment of anharmonic effects. The bonding 
of anionic hydrogen-bonded complexes is discussed and contrasted with that of 
neutral hydrogen-bonded complexes. 


1 NATO/NSF Postdoctoral Fellow. Mailing Address: NASA Ames Research 
Center, Moffett Field, California 94035 
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Introduction 


Over the past sixty years, hydrogen-bonded complexes have attracted consid- 
erable attention from chemists. Much of the interest has been directed at the under- 
standing of the nature of the relatively weak bonding present in neutral hydrogen- 
bonded complexes. To this end, severed different hydrogen bonding decomposition 
schemes have been developed. The basis for the classical description of hydrogen 
bonding was presented in a review 1 by Coulson in 195 1. The classical hydrogen 
bond energy is decomposed into four distinct components - (1) the electrostatic 
energy; (2) the delocalization energy (commonly referred to as the charge transfer 
energy); (3) the repulsive energy and (4) the dispersion energy. Since Coulson lim- 
ited his review to hydrogen-bonded complexes involving a polar molecule containing 
an electronegative atom (such as N, 0 or F) and a molecule containing a polar A-H 
bond (where A = N, 0 or F), the electrostatic interaction is viewed as the dominant 
attractive force. For some Van der Waals complexes Morokuma and coworkers have 
demonstrated 2-4 that the electrostatic energy may be very small or even represent 
a repulsive force. However, for most hydrogen-bonded complexes the electrostatic 
interaction will be attractive. The dispersion energy also represents an attractive 
force and thus, in Coulson's review, the “repulsive force” is the only interaction 
which separates monomers A and B. The explanation of the physical nature of 
this “repulsive force” is based, not surprisingly, on electron-electron repulsion, i.e. 
the mutual repulsion of the electron cloud of monomers A and B, and quantum 
mechanical effects are not discussed. 

Subsequently Morokuma and coworkers 2-4 extended and adapted this decom- 
position scheme into a rigorous quantum mechanical approach as viewed through 
the self-consistent-field (SCF) ab initio method. There are six components in this 
decomposition scheme - 1) electrostatic; 2) polarization; 3) exchange repulsion; 4) 
charge transfer; 5) “MIX”; and 6) “CORR.” The CORR term is the contribution 
of electron correlation which Morokuma and coworkers did not investigate in detail 
though they stated that the most significant portion of the intermolecular correla- 
tion energy is known as the dispersion energy which is an instantaneous effect due 
to the simultaneous correlation of electrons in monomer A and monomer B. The 
MIX term is the higher order couplings of the first four components. The polariza- 
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tion interaction is the distortion of the electron density of A (B) by the presence 
of monomer B (A) and higher order effects. In applying this decomposition scheme 
to normal hydrogen-bonded complexes, Umeyama and Morokuma 4 concluded that 
the binding in these chemical systems is mostly electrostatic in nature with a small 
but significant contribution from the charge transfer energy. 

More recently qualitative approaches based upon electrostatic and polariza- 
tion interactions have been developed for the theoretical prediction of molecu- 
lar structures 5,6 and vibrational frequency shifts 7 of hydrogen- bonded complexes. 
When applied to neutral hydrogen-bonded complexes, both of these methods yield 
qualitatively correct results, although their accuracy is generally not quantitative 
and, in some cases, not even semi-quantitative 5 " 7 . Furthermore, based upon the 
results of these approaches which have been reported thus far, it seems likely that 
these classical, perturbative approaches will break down as the binding energy of 
the complex increases. Since anionic hydrogen-bonded complexes are typically much 
more strongly bound, these simple approaches are not likely to be as successful. 

A more rigorous approach to the study of weakly bound systems which has 
been applied with much success has been the use of ab initio quantum mechanical 
methods 3 . Numerous studies have demonstrated that the SCF method (coupled 
with a large one-particle basis set) is capable of describing a weak hydrogen bond 
reasonably well, except for the dispersion energy. However, the recent formulation 
and development of better and more efficient electronic structure methods has en- 
abled the direct quantum mechanical investigation of weakly bound molecular com- 
plexes at correlated levels of theory. For example, Handy and coworkers 9 " 1 - have de- 
termined the equilibrium structures of several weakly bound hydrogen-bonded com- 
plexes [HCN-HF, HCN-HC1, (C 2 H 2 ) 2 , (C 2 H 2 ) 3 , FH-CO and FH ' NNO] using 
large one-particle basis sets in conjunction with second-order M 0 ller-Plesset pertur- 
bation theory (MP2) and have found good agreement with experiment. However, as 
Rice, Lee and Handy (RLH) have demonstrated 12 with their study of H 2 C0 "HC1, 
MP2 is not always adequate, especially when electron correlation effects are very 
important in the binding of the complex. RLH found that the theoretically more 
complete coupled-pair functional (CPF) approach gives much better H 2 CO and HC1 
monomer properties and, consequently, the H 2 CO" HC1 structure is in excellent 
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agreement with the limited experimental data. In particular the dipole moment 
of HC1 is much better described with the CPF approach, supporting the thesis 
that electrostatic interactions are important in hydrogen bonding. Nonetheless, the 
substantial differences in the equilibrium geometry of the H 2 CO " HC1 complex ob- 
tained at the SCF, MP2 and CPF levels of theory demonstrate the importance of 
the dispersion energy. 

Theoretical and experimental studies of anionic hydrogen-bonded complexes 
are more recent, especially in the gas phase. With the aid of three theoretical 
studies 13-15 , Kawaguchi and Hirota 16 have recently detected and analyzed the first 
high resolution infrared (IR) band of an anionic hydrogen-bonded complex (FHF - ). 
There have been several theoretical studies of anionic hydrogen-bonded complexes, 
though very few of these have determined equilibrium structures and molecular 
properties beyond the SCF level of theory. Furthermore, none of the theoretical 
investigations have studied the decomposition of the hydrogen bond energy of an 
asymmetric anionic hydrogen-bonded complex. Umevama et al. have performed 1, 
a decomposition of the hydrogen bond energy of FHF - and find, not surprisingly, 
that charge transfer is much more important than for neutral hydrogen-bonded 
complexes. However, the decomposition analysis of FHF - is almost certainly not 
representative of asymmetric systems since FHF - adopts a equilibrium struc- 

ture. 

In some respects anionic hydrogen-bonded complexes provide more of a chal- 
lenge than neutral and cationic hydrogen-bonded complexes for both experimental- 
ists and theoreticians. For example, the high resolution IR spectroscopist must deal 
with the very small population of anions that can be generated. Moreover, once 
a sufficient population has been attained, the analysis of the spectrum is further 
complicated by the presence of many other ionic species. The difficulty in the ab 
initio study of anionic species is well documented (see for example references 18-23). 
This difficulty generally arises due to the greater importance of electron correlation 
in anionic species. 

However, in other respects the study of anionic complexes is much easier than 
the study of similar cationic complexes. From an experimental viewpoint, the large 
binding energies of anionic hydrogen-bonded complexes should make their gener- 
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ation an easier task. In order to assess the implications of theoretical studies of 
anionic complexes, consider H 5 0^ and H3O2 • These two systems are isoelectronic 
and it is likely that electron correlation effects will be more important for H3O2 
than for HsO^. However, H3O7 has two fewer nuclei and therefore six fewer nu- 
clear degrees of freedom. The significant point is that most of the nuclear degrees 
of freedom which have been eliminated are large amplitude motions, an adequate 
treatment of which requires knowledge of a large portion of the potential energy 
surface (PES) as well as a sophisticated treatment of the nuclear motion problem. 
Existing methods for the accurate determination of the vibrational energy levels of 
polyatomic species which go beyond the harmonic oscillator approximation and are 
capable of adequately treating large amplitude motions are highly dependent upon 
the number of large amplitude nuclear degrees of freedom. Thus, while the de- 
scription of the electronic structure of anionic hydrogen-bonded complexes is more 
difficult, the accurate solution of the nuclear motion problem should be more feasi- 
ble. Therefore, the results of the current study provide data which will ultimately 
enable the detailed theoretical investigation of molecular systems with several large 
amplitude nuclear degrees of freedom. 

For anionic hydrogen-bonded complexes, the most difficult region of the PES 
to describe theoretically is the proton transfer coordinate which corresponds to the 
process AH -f B" — ► A - 4- HB. The difficulty arises due to the possible exis- 
tence of two minima corresponding to A~"‘HB and AH"'B . Several studies have 
investigated this region of the PES for symmetric and asymmetric anionic hydrogen- 
bonded complexes (see for example references 24-27). In addition, one of these" 7 
also examined the adequacy of various vibrational analysis techniques. These stud- 
ies have demonstrated that electron correlation effects are vitally important 27 in 
obtaining a reliable description of the PES along the proton transfer coordinate, and 
that when A and B are both very electronegative atoms there is generally no barrier 
(and hence no second minimum) to proton transfer for asymmetric systems 24-27 . 

To date, no high level theoretical investigations of anionic hydrogen-bonded 
complexes involving HF and CN~ or H2O and CN — have been reported. Experi- 
mentally, the IR spectrum of the M + FHCN~ ion pair (M + being an alkali metal 
cation) has been studied via matrix isolation techniques by Ault 28 . Fundamental 
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vibrations were observed in the 1100 cm -1 , 1800 cm -1 and 2500 cm -1 regions 
and assigned to a bending mode, the proton transfer mode (mostly H-F stretch) 
and the C-N stretching mode, respectively. Ault also observed that the 1100 cm -1 
band split into two components which were attributed to the presence of the metal 
cation M + . All three modes varied somewhat depending upon the composition 
of the matrix and the reactants used to form the M + FHCN - ion pair. We note 
that Ault does not seem to have considered the existence of the FHNC - isomer. 
Larson, McMahon and Szulejko 24,29 have determined the binding energies (i.e. hy- 
drogen bond strength) of both the FH 'CN - and the HjO'CN - complexes. The 
former is more strongly bound (21.1 kcal/mole) while the latter’s binding energy 
(12.7 kcal/mole) is still much larger than that of a typical neutral hydrogen-bonded 
complex. 

The purpose of this study is to obtain a better understanding of the electron 
correlation requirements in the ab initio study of asymmetric anionic hydrogen- 
bonded complexes and a more complete understanding of the nature of the bonding 
present in such systems. Thus, the conclusions of the present study will be useful 
in deciding upon the level of ab initio theory necessary to determine accurately the 
PES of an anionic hydrogen-bonded complex. The theoretical approach is described 
in the next section. The following sections contain a presentation and evaluation of 
our results. Concluding remarks are presented in the final section. 


Theoretical Approach 

It is well known 18-21 that large basis sets are necessary in order to obtain 
highly accurate results for anionic systems. Therefore, a single, large one-particle 
basis set has been used in this study. This basis consists of Dunning’s 30 (5s3p) 
contraction of Huzinaga’s 31 [10s6p] Gaussian primitive set for the heavy atoms (C, 

N, O and F). For hydrogen, the standard (3s) contraction 30 of the [5s] primitive 
set 31 was used. The hydrogen s function exponents were scaled by a factor of 1.49, 
as suggested by Dunning. In order to describe better the anionic nature of these 
systems diffuse s and p functions were added to the heavy atom basis (a J)P (C) = 

O. 04812, 0.03389; a,, p (N) = 0.06742, 0.04959; a,, p (0) = 0.08993, 0.05840; a JiP (F) = 
0.1164, 0.07161) while a diffuse s function was included in the hydrogen atomic basis 
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(a,(H) = 0.06696). These orbital exponents were determined in an even tempered 
manner using a method suggested previously 18 . Finally, two sets of polarization 
functions were added to all the atomic basis sets. The orbital exponents of the d 
polarization functions are ad = 1.5, 0.35 for the heavy atoms and a p = 1.4, 0.25 for 
hydrogen. These are the values suggested by van Duijneveldt 32 and have been used 
previously 18 in the study of anionic systems. This basis set is designated TZ2P -r 
diffuse. In all cases, the full complement of six Cartesian d functions was included 
in the basis giving 110 basis functions for the two larger complexes and 100 basis 
functions for the two smaller complexes. Linear dependency tests of the one-particle 
basis set were performed routinely and no problems were encountered. 

The first ab initio method utilized is the simplest, namely the restricted 
Hartree-Fock (RHF) SCF technique. As discussed previously, electron correlation 
effects must be included in order to account for the dispersion energy in hydrogen- 
bonded complexes. Moreover, it is reasonable to expect that the dispersion energy 
(or electron correlation effects) will be more important for anionic hydrogen-bonded 
complexes because of the diffuse, polarizable nature of the electron cloud of anions. 
Therefore, three different electron correlation methods have been used in order 
to investigate the importance of electron correlation effects. The first approach 
is second-order Mpller-Plesset perturbation theory 33 (MP2). The second, another 
commonly used method, is singles and doubles configuration interaction (CISD) 
which is based upon the variational principle. The third, and theoretically most 
complete method, is the singles and doubles coupled cluster approach (CCSD). The 
MP2 and CCSD methods have the advantage of being exactly size extensive and size 
consistent 34 . The CISD technique includes configuration mixing which MP2 does 
not take into account, but CISD is an n 6 procedure (MP2 is an n 5 procedure) where 
n is the number of active molecular orbitals. The CCSD method does allow config- 
uration mixing but is somewhat more expensive than CISD 22 . Also, the energetics 
of complicated chemical reactions are more easily computed using size consistent 
methods since “super-molecule” energies are not necessary. Therefore, based upon 
the above discussion and previous results 35 the CCSD method is expected to yield 
the most reliable results. 

Equilibrium structures of the complexes have been obtained with the SCF, 
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MP2 and CISD methods. Due to the computational cost and the available compu- 
tational facilities, at the CCSD level of theory it was only possible to optimize the 
FH "CN“ and FH'"NC~ complexes. However, single point CCSD energies at the 
MP2 and CISD equilibrium structures have been performed for the H 2 0 'CN and 
H 2 0 "NC~ pair of dimers. Also, in order to reduce the CISD expansions the heavy 
atom ls-like core molecular orbitals were required to be doubly occupied in all 
configurations and the corresponding virtual counterpart was deleted from the pro- 
cedure. The same procedure was also used in the CCSD optimization of the two 
FH;CN _ complexes. 

As noted in the introduction, the structures of many neutral hydrogen-bonded 
complexes are strongly dependent upon the respective monomer properties. Thus, 
in an attempt to judge better the reliability of the theoretical predictions of the 
complexes, all possible monomers have been studied using the basis set and ab 
initio methods described above. 

In most cases, analytic energy gradient methods 36 have been employed to lo- 
cate precisely the equilibrium structures. Analytic energy second derivative meth- 
ods have been used to determine the SCF 37 and MP2 38 ’ 39 Hessian matrices while the 
CISD and CCSD Hessians were obtained numerically by taking central differences 
of analytic gradients. Infrared intensities have been determined via the double har- 
monic approximation. The dipole derivatives were determined analytically at the 
SCF level of theory and central differences of dipole moments were utilized at the 
CISD and CCSD levels of theory. In all cases, dipole moments were determined with 
respect to the center of mass and evaluated as energy derivatives 40 . In the numer- 
ical central difference procedures, energy invariance relationships for the Hessian 41 
and dipole derivative 42 matrices were used in order to reduce the number of gradi- 
ent evaluations. This is the most efficient numerical procedure for the evaluation of 
dipole derivatives provided that the numerical Hessian is also required. 

All SCF and MP2 investigations were performed with the Cambridge Analytic 
Derivatives Package 43 (CADPAC), while the CISD 44 - 45 and CCSD 46 ’ 47 studies were 
performed with the Berkeley suite of programs modified to run on a Cray X-MP. 
The CCSD studies of the dimers were performed with a recently developed 22 vector- 
ized CCSD method. SCF, MP2, CISD and monomer CCSD calculations were per- 
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formed at the University of Cambridge. The CCSD optimizations of FH"*CN and 
FH’"NC“ were performed on the Cray X-MP/48 at NASA Ames Research Center. 


Monomer Properties 

The equilibrium structures, total energies and dipole moments of the various 
monomer fragments are presented in Table 1. The CCSD method provides much 
better agreement with experimental structures and dipole moments. In fact, for 
the neutral molecules the magnitude of the errors in the CCSD prediction of the 
equilibrium structures are less than half those present for the CISD and MP2 struc- 
tures. The only exception arises for the bond angle in H 2 O where the MP2 result is 
fortuitously in better agreement with the experimental value. However, the CCSD 
value is only 0.3° too large. Nonetheless, the main conclusions to be drawn from 
the results of Table 1 are that the CCSD method, as expected, performs better 
than either the CISD or MP2 approaches and, more importantly, that for many 
chemical systems quantitatively accurate structures (i.e., Ar c < 0.001 Aand A0 e 
< 0.5° ) may be obtained with the CCSD electron correlation procedure coupled 
with a large one-particle basis set. 

Another important aspect concerns the equilibrium structure of OH*. Note 
that the absolute magnitude of the error of the CCSD bond length is significantly 
larger than for the A-H (A=C,O.N) bonds of the neutral molecules. The electron 
correlation energy for anions is generally larger than for isoelectronic neutrals (as 
evidenced here by comparing the correlation energies of OH” and H 2 O). Thus, 
anions usually require a more rigorous treatment of electron correlation in order to 
obtain accuracy comparable to that obtained with neutral molecules. With these 
considerations, it is not too surprising that the CCSD bond length of OH” is not 
as accurate (compared to experiment) as the CCSD O-H bond length in H 2 O. 

Based upon previous experience, the only other geometrical parameter which 
is potentially difficult for ab initio methods is the C-N triple bond present in CN”, 
HCN and HNC. The results given in Table 1 confirm the inherent difficulty in 
adequately treating the C-N triple bond, though again the CCSD equilibrium values 
for HCN and HNC are superior to either the CISD or MP2 quantities. For these 
three molecules, the C-N bond distance decreases in the order CN“ > HNC > HCN, 
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which yields insight into the nature of the carbon and nitrogen lone pairs. Since the 
C-N bond distance in CN - is longer than that in C-N radical, the HOMO exhibits 
slightly anti-bonding character. Thus, when a proton is attached to form either 
HCN or HNC the C-N bond distance shrinks due to the polarization of electron 
density away from the C-N linkage. Therefore, since the C-N distance in HCN is 
shorter than that found in HNC, we may conclude that C contributes more to the 
antibonding characteristics than does N. As an aside to the above discussion it is 
interesting to note that the CISD and CCSD correlation energies for these three 
molecules increase in magnitude in the order HNC < CN - < HCN. At the MP2 
level of theory the correlation energy of CN - is slightly larger than that for HCN. 

The dipole moments of the monomers are predicted to almost equal accuracy 
with the CISD or CCSD methods, though the CCSD dipole moment is usually 
in better agreement with experiment. A noteworthy point which has particular 
relevance to this study is the fact that the rather sizable dipole moment of CN - 
(0.64 D) has carbon at the negative end. Consequently, there are competing effects 
as to which lone pair of electrons (the C lone pair or the N lone pair) will act as 
the better Lewis base, or in other words, which end of CN - will form the stronger 
hydrogen bond? The nature of these phenomena may be understood by considering 
electron density maps of the C and N lone pair molecular orbitals which have been 
given by Taylor et al. i7 The C lone pair orbital is broad and diffuse whereas the 
N lone pair orbital is tighter. Thus, the C lone pair electrons will produce a larger 
attraction on the proton of the hydrogen-containing monomer in the AH" CN - 
complex whereas for the AH ' NC - complex the N nucleus will have a stronger 
interaction with the electron cloud of the hydrogen-containing monomer since it 
will be able to approach more closely (i.e. form a shorter hydrogen bond). Thus, it 
is not evident, a priori, which of the two complexes will be more stable. Therefore, 
if CN - is one of the monomers of a hydrogen-bonded dimer, it will be necessary to 
investigate both AH ' CN - and AH ' NC - . As we shall demonstrate, for AH being 
either HF or H 2 O both sets of isomers are nearly isoenergetic. 

The harmonic vibrational frequencies and infrared (IR) intensities of the 
monomers are reported in Table 2. As has been noted by several authors recently, 
the CISD method seems incapable of properly describing the curvature of the PES 
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around an equilibrium point whereas the size extensive MP2 and CCSD methods 
both yield quite good harmonic vibrational frequencies. The CISD harmonic fre- 
quencies are consistently too high even with the rather large one-particle basis set 
used in this study. This particular inadequacy with CISD is believed to be related 
to the lack of size extensivity 49 , though no direct proof has as yet been given. 

Somewhat surprisingly the MP2 and CCSD harmonic frequencies are about 
equally accurate for this set of molecules (with the TZ2P+diffuse basis set), with 
the CCSD harmonic frequencies for OH - , HCN and HNC being somewhat better 
than the MP2 values and the CCSD harmonic frequencies for HF and H 2 0 being 
marginally worse than the MP2 quantities. In any case, the MP2 and CCSD har- 
monic frequencies for the monomers are in very good agreement with experimental 
values with the possible exception of the bending mode in HCN and HNC. This 
particular normal mode is very sensitive to specific basis set deficiencies and the 
interested reader is referred to references 35 and 50 for more details of this effect. 
Based upon the CCSD and experimental results of the C-N stretching normal mode 
of HCN and HNC, the experimental harmonic frequency of CN - can be estimated 
to lie near 2076 cm -1 . Using the u e x t = 11.3 cm -1 determined by Taylor and 
coworkers 48 , the experimentally unknown fundamental frequency is predicted to lie 
at 2053 cm -1 . This value is in excellent agreement with the high level calculations 
of Botschwina 51 (2052 ± 6 cm -1 ). 

The C-N stretch harmonic frequency decreases in the order HCN > CN 
> HNC. Based upon the previously discussed C-N bond distances, the CN har- 
monic frequency would probably have been expected to be the lowest. This result 
demonstrates that caution must be exercised in relating geometric and vibrational 
properties. 

The IR intensities reported in Table 2 are consistent 52 with the expectation 
that electron correlation tends to reduce the magnitudes. The CCSD IR intensi- 
ties demonstrate that while CISD IR intensities are a vast improvement over SCF 
quantities, the CISD procedure still underestimates the correlation contribution to 
IR intensities. This observation is entirely consistent with a recent study 49 on the 
effects of triple and quadruple excitations in the Cl electron correlation procedure 
where it was shown that, like the electronic energy, many molecular properties 
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tend to converge from one direction as the excitation level is increased (i.e., do not 
exhibit oscillatory convergence). Thus, based upon the IR intensity and dipole mo- 
ment data, we may conclude that the CCSD approach better describes the electrical 
properties of the molecular systems included in this study. 


Energetics 

In order to determine which dimer will represent the global energy minimum 
(e.g. FH-CN" or F""HCN), it is necessary to consider the enthalpy of the two 
reactions 

HF + CN" + AH 298 -> F" + HCN (1), 

and 

H 2 0 + CN" 4- AH 2 98 OH" + HCN (2). 

If AH 298 is positive, then the dimer will correspond to the reactants. Table 3 
contains ab initio and experimentally derived values for AH 29 B- The data in Table 
3 clearly indicate that the reactants of equations (1) and (2) should form the more 
stable dimer. This situation arises due to the large electron affinity of CN and the 
large F-H and 0-H bond energies. Therefore, as discussed above, we must consider 
two sets of isomers corresponding to hydrogen bond formation through the C or N 
end of CN“. 

Table 4 contains the total energy, binding energy and dipole moment of each 
complex, determined at the equilibrium structure for the given level of theory. The 
most important point to notice is that both sets of isomers are nearly isoenergetic. 
Thus, it is not possible to say definitively which stationary point represents the 
lowest energy structure. However, the ab initio data are reliable enough to con- 
clude that the actual difference between the isomer’s binding energies will not be 
greater than 5 kcal/mole. Thus, depending upon how the complexes are formed, 
it is possible that both isomers will be present under a given set of experimental 
conditions. 

The binding energies in Table 4 have incorporated a correction for the basis 
set superposition error (BSSE). The BSSE was determined using the counterpoise 53 
method at the SCF and CCSD levels of theory . It is well established 8 that the BSSE 
is generally larger at a correlated level of theory and that in order to reduce the 
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BSSE at a correlated level of theory a very large one-particle basis set is required 8 ’ 10 
. Therefore, we include the CCSD BSSE for the energetics determined with electron 
correlation methods. 

The best theoretical estimates of the [^O'CN - ; H20"'NC _ ] binding en- 
ergy (14.5 and 14.7 kcal/mole, respectively) are in good agreement with the 
experimental 29 value, 12.7 ± 0.8 kcal/mole. The agreement between theory and 
experiment for the binding energy of the [FH"’CN - ; FH'NC - ] set of isomers is 
also good, again being somewhat too large. However, the difference between the- 
ory and experiment is somewhat larger for the FH;CN~ pair of complexes. The 
experimental value may be somewhat too low for this complex and support for this 
assertion is found by comparing the theoretical and experimental binding energies 
of the F~ 'H 2 0 complex. A similar 27 level of theory to that used in this study gave 
a binding energy of 23.2 kcal/mole for F - "H 2 O with the experimental quantity 
being 23.3 kcal/mole. Thus, the results of this study suggest that the FH;CN~ pair 
of complexes may be slightly more strongly bound than F - '"H 20 , but experimental 
values suggest the opposite situation. As we shall show, other molecular properties 
determined via ab initio methods (such as the IR intensity of the proton trans- 
fer mode) are consistent with the FH;CN~ pair of complexes being more strongly 
bound. 

For both sets of isomers the complex which is hydrogen bonded through the 
N end of CN~ has a much larger dipole moment. This situation occurs because, in 
all cases, the negative end of the dipole moment of the dimer is the CN~ end of the 
complex and so the dipole moment of the complex is greater in magnitude when the 
negative end of the CN - moiety is farthest from the center of mass of the dimer. In 
addition, the difference between the dipole moment of FH’NC - and FH ’CN - is 
0.63 D (CISD), almost exactly the CISD dipole moment of CN - (0.60 D). 

As is usual for hydrogen-bonded complexes, the dipole moment of the complex 
is greater than the vectorial sum of the two monomers. However, in this case the 
increase is much larger than normal and is partially due to the charged nature of the 
complex and the large change in the relationship between the center-of-mass and 
the center-of-electron charge which occurs upon formation of the complex. Also, 
the large polarizability of CN - probably contributes to the sizable dipole moment 
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of the complexes due to polarization of the CN electron cloud away from the HF 
or H 2 O species. 

Structures 

Table 5 lists the ab initio equilibrium structures of the FH'CN" and 
FH "NC" anionic hydrogen-bonded complexes, with those of the H20"'CN and 
H 2 0"NC" dimers presented in Table 6. See Figures 1 and 2 for the definition of 
the geometrical parameters contained in Table 6. 

For the FH"CN~; FH "NC" pair of complexes the MP2 level of theory greatly 
overestimates the effects of electron correlation and CISD underestimates the impor- 
tance of electron correlation; consistent with the results obtained for the monomers. 
The CCSD method predicts bond lengths which are between the MP2 and CISD 
values, but which are much closer to the CISD values than the MP2 quantities. 
This indicates the importance of electron correlation. 

Comparing the complex tfh and r cn with the monomer bond lengths we note 
that the H-F bond distance increases, as is typical upon hydrogen bond formation, 
but that the C-N bond distance decreases relative to CN . This effect may be due 
to the loss of some of the C-N antibonding character and is supported by the earlier 
observation that the C-N distance in CN" is longer than in either HCN or HNC. 
This explanation is also consistent with the experimentally observed blue shift in 
the C-N stretch frequency of HCN’HF 54 . Since electron density is drawn away 
from the C-N bond, the C-N antibonding character is reduced leading to a shorter 
C-N distance in the complex. The C-N stretch frequency of HCN often exhibits a 
blue shift in neutral hydrogen bonded complexes 55 . 

Interestingly, although the heavy-atom distance (Rf..c or Rf..n) I s smaller for 
FH' "NC" (due to the shorter hydrogen bond), the C-N distance is more affected 
(relative to CN") in FH"’CN". This result tends to suggest that C contributes 
more to the C-N antibonding characteristics. Another noteworthy feature of the 
heavy-atom distances, is that for both isomers the CCSD level of theory predicts 
the largest distances while the MP2 method dramatically underestimates the heavy- 
atom lengths. 

For the H 2 0"CN" and H 2 0"NC" pair of complexes, CCSD geometry opti- 
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mizations were not possible. However, based upon the above comparisons between 
CCSD, CISD and MP2 for the FH;CN - pair, it is reasonable to expect that the 
CCSD equilibrium structures will be intermediate between the CISD and MP2 opti- 
mum geometries and probably somewhat closer to the CISD structures. The various 
equilibrium structures of the H 2 0"CN - and H 2 0"NC - anionic complexes given 
in Table 6 exhibit tendencies similar to those reported above for the FH;CN - pair 
of dimers. The O-Hj bond distance (where Hj is involved in the hydrogen bond, 
see Figures 1 and 2) elongates upon complexation while the C-N linkage decreases 
relative to that in CN - . In addition, the N-Hi hydrogen bond distance is again 
shorter than the C-Hi hydrogen-bond distance; in this case by 0.138 A. The differ- 
ential heavy atom distance (i.e., Rco - Rjvo = 0.1402 A, CISD) is also larger than 
for the FH;CN - pair. 

Unique to the H 2 0;CN - complexes is the decrease in the 0-H 2 bond distance, 
the closing of angle 7 (see figures 1 and 2 for the definition of 7) and the non-linear 
hydrogen bond (i.e., A-H 'B do not lie in a straight line). The decrease in the 0-H 2 
bond distance seems natural due to the longer O-Hi distance, though, this result 
seems to imply that the electron density of the H 2 0 monomer unit is polarized 
towards the CN - monomer unit. While this phenomenon would be expected for 
neutral hydrogen-bonded complexes, it is not necessarily expected for the case where 
one of the monomers is an anion. However, the shorter 0-H 2 distance may be 
related to the decrease in the O-H-O angle 7. In other words, long range attractive 
forces between the H and the electron cloud around the C and/or N will result 
in a decrease in both r o-h 2 and 7. Such long range attractions also explain the 
non-linear hydrogen bond. 

Aside from the hydrogen bond distance and the associated heavy atom dis- 
tance, the main structural difference between H 2 0"CN - and H 2 0"NC - is the 
angle a. The smaller angle a for H 2 0 "NC - represents a larger deviation from 
linearity and is consistent with long range attractive forces between H 2 and the 
electron cloud around C (N in the case of H 2 0'"CN - ). Since the electron density 
around the C end of CN - is more diffuse, there is a stronger interaction between 
H 2 and C in H 2 0 "NC - than between H 2 and N in H 2 0"CN - . Thus, the angle a 
is smaller by about 4° for H 2 0'”NC - . 
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Diagonalization of the mass-weighted Hessian matrices explicitly demonstrate 
that each of the four complexes represents a true minimum on the PES. Therefore, 
attempts were made at the SCF level of theory to locate the transition structure 
between FH 'CN - and FH"NC~. However, due to the nature of interactions be- 
tween two closed-shell monomers the potential energy surface is very flat in this 
region in several degreees of freedom, though the total energy does rise as the CN - 
moiety rotates. Performing the full geometry optimization is somewhat complicated 
and so the actual stationary point structure of the transition state was not pursued 
further. 


A search of the potential energy surface along the proton transfer coordinate 
was also performed in order to determine whether a second minimum (corresponding 
to F“"HCN or P~-HNC) exist. The search along the PES in this coordinate is 
significantly easier since all the atoms were constrained to be collinear. However, a 
second minimum (and corresponding transition state) could not be located. A brief 
discussion of the nature of the PES along the proton transfer coordinate is in order. 
Generally a transition state (TS) on a PES arises due to an avoided crossing of 
two states of the same symmetry. Thus, the SCF method is often not an adequate 
reference function for the TS. However, there are many different types of avoided 
crossings and in this particular case the SCF wavefunction should be a reasonable 
reference. This situation arises due to the fact that the orbital occupations of the 
reactants (A _ ’"HB) and products (AH 'B - ) are the same. What does occur as 
the proton is transferred is that two reactant molecular orbitals (a lone pair MO 
on A~ and a bonding MO in HB) change character and become two product MO’s 
(a lone pair on B - and a bonding MO in AH). However, since these MO’s are 
of the same symmetry, the transition from reactant to product MO’s is smooth 
along the proton transfer coordinate. We note that the reactant and product MO’s 
belonging to the same irreducible representation is a necessary but not sufficient 
condition for a smooth transition. Nevertheless, in this specific case the transition 
from reactant to product MO’s appears to be smooth. Therefore, the SCF function 
should represent a reasonable reference from which to evaluate dynamical electron 
correlation effects. 
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Vibrational Spectra 

The harmonic vibrational frequencies and infrared intensities for the four an- 
ionic hydrogen-bonded dimers included in this study are presented in Table 7. The 
experimental fundamentals which Ault 28 measured in matrix isolation IR studies 
are included for comparison, though, because of the rather large anharmonicities 
which the stretch modes are expected to exhibit, near quantitative accuracy with 
harmonic frequencies is not possible. The most astonishing result from Table 7 is 
the variation of the harmonic frequencies u,, FH~CN“ and u, a FH-NC~ with re- 
spect to level of theory. Note that the normal mode associated with u> a corresponds 
to the proton transfer coordinate (i.e., AH 'B — ♦ A 'HB), which for the FH,CN 
complexes is predominately the H-F stretch. Quite clearly, an adequate treatment 
of electron correlation is extremely important in properly describing the shape of 
the potential energy surface along this coordinate. Interestingly, the large varia- 
tions in (e.g., for FH "CN“, 3320 cm' 1 SCF, 2352 cm" 1 MP2 and 2844 cm' 1 
CISD) would probably not have been predicted based upon the different equilib- 
rium H-F bond distances (0.947 A SCF, 1.010 A MP2, and 0.977 A CISD), though, 
not surprisingly, there is a strong correlation between the H-F bond distance and 
the harmonic frequency. In fact, the nearly linear relationship between Wjand Tfjjr 
allows the CCSD oij to be estimated as ~ 2689 cm -1 for FH'CN and ~ 2912 
cm" 1 for FH-NC-. 

Given the large variation of with level of theory, it may seem very difficult 
to arrive at a reliable theoretical prediction for the fundamental band center V\. 
However, studies on similar systems have demonstrated that the individual har- 
monic frequency and anharmonic correction quantities converge much more slowly 
(with respect to level of theory) than does the combination, i.e., the fundamental 
band center. For example, in the study 13 of FHF - by Janssen et al. the harmonic 
frequency of the antisymmetric stretch u >3 varies from 627 cm -1 to 1538 cm -1 while 
the fundamental u 3 varies only from 1427 cm -1 to 1703 cm -1 . A possible expla- 
nation for this observation may be that the A"'H"'B system should be viewed as a 
particle in a one- dimensioned box, where the distance Rab defines the box in which 
the proton is allowed to move. Thus, we may expect the distance Rab to converge 
more quickly (with respect to level of theory) than Tab and tbh- ! n reexamining 
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the theoretical structures in Table 5 we note that once an iterative electron correla- 
tion procedure is used, then the above conditions are met (i.e., AR < Ar). In any 
case, the above explanation seems feasible and will no doubt be scrutinized as more 
theoretical studies concerned with the prediction of the fundamental vibrational 
frequencies of this type of system are performed. 

The second most striking feature of the IR spectrum of the FH;CN~ pair 
of complexes is the extremely large intensity exhibited by <*>i. Though a large 
IR intensity is expected for a mode which corresponds to proton transfer, the IR 
intensities of Wj FH "CN" and uq FH "NC _ are even larger than the IR intensity 
reported for the analogous mode of F - 'H 2 0. However, the IR intensity reported 
for the asymmetric stretch of FHF~ is substantially larger than the u>i FH;CN“ 
quantities. Interestingly, there appears to be a direct correlation between the IR 
intensity of the proton transer mode and the binding energy of the dimer. The 
appropriate IR intensity and the ab initio binding energy both decrease in the 
order FHF" > FH;CN" > F - " H 2 0 > H 2 0;CN _ . This correlation suggest that 
the larger the anionic dimer binding energy then the flatter the potential energy 
surface along the proton transfer coordinate leading to a larger amplitude motion. 
Of the remaining FH;CN _ normal modes, W 4 and possibly W 3 should be observable 
with u > 2 of FH-NC" also a possibility. The IR intensities of the H 2 0;CN~ pair of 
complexes are more evenly distributed and, therefore, there are several vibrational 
modes which should be observable. 

For the H 2 0;CN - pair of complexes the proton transfer vibrational mode is 
u> 2 - The variation of u ; 2 with respect to level of theory is much smaller than was 
exhibited by the FH"’CN - and FH'"NC~ pair, though it is still substantial. For 
example, u ; 2 for H 2 0-CN- is 3786 cm " 1 , 3174 cm " 1 , and 3497 cm " 1 for the SCF, 
MP2 and CISD levels of theory, respectively. The smaller variation of u > 2 H 2 0;CN - 
relative to u>j FH;CN - was, however, to be expected due to the smaller binding 
energy of the H 2 0;CN~ pair. 

Another manifestation of the smaller binding energy of the H 2 0"' CN~, 
H 2 0"NC~ pair of complexes is the lower C-N stretch harmonic frequency rela- 
tive to the FH ' CN - , FH "NC - pair. As noted previously, the harmonic frequency 
of CN~ is less than u > 2 in HCN. Thus, in an analogous manner the lower C-N stretch 
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frequency in the H 2 0;CN~ pair is consistent with a smaller interaction between the 
H 2 0 and CN~ monomers than exists between the HF and CN monomers. These 
observations also suggest that there is a smaller degree of charge transfer in the 
H 2 0;CN~ complexes than present in the FH;CN" pair. However, the above obser- 
vations do not indicate the relative importance of charge transfer in the bonding 
mechanism. 

Not surprisingly, the C-N stretch harmonic frequency of all four complexes 
exhibits a blue shift relative to the C-N stretch in HCN, CN" and HNC. As discussed 
earlier, neutral HCN hydrogen-bonded complexes (such as HF'HCN) often exhibit 
a blue shift in the C-N stretch due to the loss of C-N antibonding character upon 
complexation. 

By comparing the vibrational spectra of the two FH;CN complexes or the 
H 2 0;CN - pair, it is evident that it would be difficult to distinguish between the 
two isomers based upon the vibrational frequencies alone. However, due to the 
differences in their structures, the best method of distinguishing the two isomers 
will be via analysis of a ro- vibrational band. The rotational constants presented in 
Table 8 confirm this hypothesis since the differences are well within the accuracy of 
high resolution spectroscopy. 

In order to predict accurately the fundamental band centers of the vibra- 
tional modes of these complexes, a potential energy function including very high 
orders (e.g., octic terms) in some of the degrees of freedom would be required. In 
addition, a high level approach to the solution of the nuclear Schrodinger equa- 
tion, which explicitly accounts for large anharmonic couplings, would be necessary. 
This procedure would obviously be very expensive and is beyond the scope of the 
present study. However, the vibrational analysis that we have performed has lead 
to further insight concerning the proton transfer vibrational mode, the most likely 
fundamental of the FH;CN — and H 2 0;CN - anionic dimers to be experimentally 
observed. Furthermore, the similarity of the vibrational spectra of the pairs of iso- 
mers has been explicitly demonstrated and a method by which the isomers may be 
spectroscopically distinguished has been noted. 


Bonding 
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The much larger binding energies found in anionic hydrogen-bonded com- 
plexes (relative to neutral hydrogen-bonded complexes) lead to questions concerning 
the nature of this interaction. For example, if Morokuma and coworkers ’ hydrogen 
bond energy decomposition scheme is applied, which components exhibit signifi- 
cantly different characteristics for anionic complexes? As discussed earlier, such an 
analysis has been performed 17 on FHF - , however, it seems likely that an asym- 
metric anionic hydrogen-bonded complex will possess quite different characteristics 
than FHF - where charge transfer is clearly very important. Moreover, the binding 
energy of FHF - (~ 39 kcal/mole) is significantly larger than that for the complexes 
included in this study. 

The three hydrogen bond components which one might intuitively expect to 
yield large attractive energies are the electrostatic, polarization and charge trans- 
fer interactions. ■ We will not discuss the electrostatic interactions here except 
to note that the detailed structure of this interaction must be very different for 
FH-CN - (H 2 0 'CN - ) and FH 'NC - (H 2 0 "NC - ) because of the reversal of the 
dipole moment of CN - . The total binding energies are very similar, however. Even 
though the dipole of CN - has been reversed, this does not mean that the total 
electrostatic energies of the two isomers are different, though it does seem probable 
that there will be a detectable difference. In the latter case, some other hydrogen 
bond energy component must compensate. 

The polarization interaction for the anionic hydrogen-bonded complexes in- 
cluded in this study must be significantly larger than exists in a similar neutral 
hydrogen-bonded complex. This conclusion is based upon the much larger polar- 
izability which anions possess (e.g., at the SCF TZ2P+diffuse level of theory the 
mean polarizability for HF, H 2 0 and CN - is 0.66, 1.14 and 3.48 A 3 , respectively). 
Furthermore, the decomposition analyses which have been performed on neutral 4 
and anionic 17 (FHF - ) hydrogen-bonded complexes provide additional support for 
this inference. 

It is difficult to assess the degree of charge transfer. One method would be 
to perform a Mulliken population analysis on the complex, and from these data 
determine the number of electrons associated with each monomer. Performing such 
an analysis on the FH 'CN - and FH“'NC - anionic complexes and comparing to 
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a similar analysis on the HCN'HF hydrogen bonded complex shows that indeed 
there is more charge transfer in the anionic species. However, as is well known, a 
Mulliken population analysis associates electrons to a given nucleus in an ambiguous 
manner. Therefore, given the rather small differences between the neutral and 
anionic complexes the validity of the results would seem to be in question. An 
alternative method would be to perform electron density difference plots between 
the complexes and their respective monomers and compare these for the anionic 
and neutral hydrogen-bonded complexes. 

Valence electron density difference plots from CISD natural orbitals have been 
performed for the FH"CN - , FH'NC - and HCN'HF hydrogen-bonded complexes 
and are presented in figures 3-5, respectively. The contour interval for all three 
plots is the same. The HCN'HF equilibrium geometry was taken from reference 
9, but the TZ2P'+ diffuse basis set of the current study was used. Noting that 
short dashed lines indicate electron depletion and solid lines indicate an increase 
in electron density, it is clear that the anionic complexes exhibit a larger charge 
transfer from the CN - species to the HF monomer than occurs in the neutral 
complex. Moreover, this conclusion is enforced by the large buildup of electron 
density behind the F atom in the anionic complexes. Interestingly, the plots also 
show a depletion of electron density in the C-N bonding region upon complexation. 
This observation is entirely consistent with earlier statements concerning the C-N 
equilibrium bond length and harmonic frequency in HCN, HNC and CN . 

Considering the above discussion and previous results 2-4,17 , a possible scene- 
rio may be suggested. It is likely that the electrostatic, polarization and charge 
transfer energy components of an anionic complex are all larger than those for a 
similar neutral complex. Moreover, as one progesses from an asymmetric complex to 
a symmetric species (i.e., the proton half way between the heavy atoms), the charge 
transfer component will become much larger. This model, then, also explains the 
large difference between the binding energies of FHF - and FH;CN - . 

Concluding Remarks 

The FH "CN - and FH 'NC - pair of anionic hydrogen-bonded complexes have 
been shown to be nearly isoenergetic and the theoretical binding energy is in good 
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agreement with experiment. The H 2 0 -CN - and H 2 0"NC“ pair of complexes 
are also very close energetically with the best ab initio binding energy again in 
good agreement with the experimental value. The equilibrium structures of the 
isomers, however, do exhibit small differences (e.g., the N"H hydrogen bond is 
shorter than the C 'H hydrogen bond) which lead to slightly different rotational 
constants. Thus, because the harmonic IR spectra of the two pairs of isomers are 
so similar, the different rotational constants provide a means by which the isomers 
may be experimentally distinguished. It is concluded, however, that an accurate 
theoretical determination of the fundamental frequencies will require a large portion 
of the potential energy surface to be investigated using a high level of electronic 
structure theory, such as CCSD coupled with a large one-particle basis set. In 
addition, a sophisticated solution of the nuclear motion problem capable of treating 
large anharmonicities will be necessary. 

Another significant outcome of this study involves the CCSD investigations of 
the monomers. This is the first study which has fully optimized molecular structures 
and evaluated several equilibrium molecular properties at the CCSD level of theory 
with a large one-particle basis set (i.e., larger than double zeta plus polarization) 
for chemical systems exhibiting a range of bonding characteristics. The CCSD 
equilibrium structures, harmonic frequencies, dipole moments and IR intensities for 
HF and H 2 0 clearly demonstrate that near quantitative results may be obtained 
for systems which are well described by a single determinant reference function. 
Although the CCSD results for HCN, HNC, and OH - have slightly larger errors, 
they are still very good and are superior to the analogous MP2 and CISD quantities. 
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Figure Captions 

Figure 1. Definition of the internal coordinates for the H 2 0-CN~ anionic hydrogen- 
bonded complex. 

Figure 2. Definition of the internal coordinates for the H 2 0"'NC anionic hydrogen- 
bonded complex. 

Figure 3. Valence electron density difference plot for the FH‘"CN~ anionic complex. 
Short dashed lines indicate a depletion of electron density while solid lines indicate 
an increase of electron density. 

Figure 4. Valence electron density difference plot for the FH‘"NC anionic complex. 
Short dashed lines indicate a depletion of electron density while solid lines indicate 
an increase of electron density. 

Figure 3. Valence electron density difference plot for the HCN'HF hydrogen- 
bonded complex. Short dashed lines indicate a depletion of electron density while 
solid lines indicate an increase of electron density. 
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Table 1 

Theoretical predictions of the total energy, optimum structure 
and dipole moment of the possible fragmentation monomers. 
Energies, bond lengths, angles and dipole moments are given in 
Hair trees, A, degrees and Debyes respectively. 


Monomer 

Structure 

Method 

Energy 

Dipole Moment 

F“ 



SCF 

-99.455226 

- 




MP2 

-99.724127 

- 




CISD 

-99.691431 

- 




CCSD 

-99.718898 

- 

HF 

r hf 

0.8985 

SCF 

-100.064599 

1.89 



0.9208 

MP2 

-100.317487 

- 



0.9149 

CISD 

-100.293761 

1.82 



0.9176 

CCSD 

-100.319608 

1.81 



0.9168 

Expt a 

- 

1.80 

OH" 

*OH 

0.9426 

SCF 

-75.413857 

1.35 



0.9655 

MP2 

-75.688586 

- 



0.9580 

CISD 

-75.654436 

1.25 



0.9624 

CCSD 

-75.685247 

1.24 



0.9643 

Expt* 

- 

- 

CN“ 

r CiV 

1.1513 

SCF 

-92.342662 

0.46 



1.1870 

MP2 

-92.695252 

- 



1.1688 

CISD 

-92.637209 

0.60 



1.1744 

CCSD 

-92.698406 

0.64 



- 

Expt 

- 

- 

h 2 o 

toh 

0.9404 

SCF 

-76.062199 

1.95 


<:hoh 

106.4° 





*oh 

0.9593 

MP2 

-76.315821 

- 


zlHOH 

104.5° 





*OH 

0.9540 

CISD 

-76.293736 

1.89 


ZHOH 

105.0° 





toh 

0.9571 

CCSD 

-76.321837 

1.88 


ZHOH 

t— * 

O 

bo 

o 





roa 

0.9578 

Expt c 

- 

1.85 


ZHOH 

104.5° 

- 




Table 1 continued 


Monomer 

Structure 

Method 

Energy 

Dipole Moment 

HCN 

*hc 

1.0572 

SCF 

-92.909715 

3.27 


r nc 

1.1235 





*HC 

1.0641 

MP2 

-93.262244 

- 


r nc 

1.1637 






1.0618 

CISD 

-93.207319 

3.08 


*NC 

1.1438 





r hc 

1.0653 

CCSD 

-93.268033 

3.03 


T NC 

1.1502 





*HC 

1.065 

Expt d 

- 

2.99 


*NC 

1.153 




HNC 

*HN 

0.9819 

SCF 

-92.892533 

2.97 


*CN 

1.1440 





*HN 

0.9959 

MP2 

-93.233324 

- 


*CN 

1.1733 





*HN 

0.9903 

CISD 

-93.185062 

3.09 


r cn 

1.1609 





*HN 

0.9939 

CCSD 

-93.244549 

3.11 


*CN 

1.1658 





*HN 

0.9940 

Expt' 

- 

3.05 


*CN 

1.1689 





a. All experimental structures refer to derived equilibrium structures. HF bond length from Ref [56] and 
dipole moment from Ref [57]. 

b. Ref [58]. 

c. Structure from Ref [59] and dipole moment from Ref [60]. 

d. Structure from Ref [61] and dipole moment from Ref [62]. 

e. Structure from Ref [63] and dipole moment from Ref [64], 




Table 2 

Harmonic vibrational frequencies and infrared intensities 
for the monomers. Frequencies are given in cm" 1 and IR intensities 
(in parentheses) in km/mol. 


Monomer 

Normal 

mode 

SCF 

MP2 

CISD 

CCSD 

Expt 

HF 

Wi(o-) 

4469 (164) 

4126 

4212 (114) 

4165 (106) 

4139 (96)“ 

CN" 

wi(cr) 

2317 (45) 

1982 

2167 (21) 

2112 (16) 

- 

OH" 

Wl(ff) 

4073 (62) 

3805 

3855 (77) 

3782 (85) 

3738 4 

H 2 0 

wi(ai) 

4130 (15) 

3841 

3919 (6) 

3865 (4) 

3832 (2)* 


w 2 (ai) 

1757 (96) 

1657 

1694 (74) 

1684 (71) 

1649 (54) 


w 3 (b 2 ) 

4233 (92) 

3967 

4021 (62) 

3972 (57) 

3942 (45) 

HCN 

Wi(ff) 

3608 (72) 

3451 

3497 (68) 

3438 (64) 

3442 (59) d 


w 2 (<t) 

2407 (11) 

2027 

2236 (2) 

2171 (0.4) 

2129 (0.2) 


w 3 (jt) 

855 (70) 

686 

734 (72) 

706 (72) 

727 (50) 

HNC 


4046 (379) 

3818 

3899 (286) 

3839 (260) 

3842* 


w 2 (<t) 

2282 (103) 

2017 

2145 (70) 

2098 (62) 

2067 


w 3 (tt) 

472 (313) 

459 

425 (277) 

442 (269) 

490 


a. All experimental frequencies are derived harmonic frequencies. The harmonic frequency is taken from ref 
[56] while the IR intensity is taken from ref [65]. 

b. The harmonic frequency is taken from ref [58]. 

c. The harmonic frequencies are taken from ref [66] while the IR intensities are taken from ref [67]. 

d. The harmonic frequencies are taken from ref [68] while the IR intensities are taken from ref [69]. 

e. The harmonic frequencies are taken from ref [63]. 


Table 3 

Thermochemical data for possible fragmentation products of the 
titled anionic hydrogen bonded complexes in kcal/mole. 


Method 

AH 298 a ’ c 

ah 298 6 

SCF 

27 

50 

MP2 

17 

36 

CISD 

21 

42 

CCSD* 

20 

41 

Expt e 

26 

45 


a) HF + CN~ + Ai?298 -*■ F~ + HCN 

b) H 2 0 + CN~ + AH 2 9s -+ 0H~ + HCN . 

c) Difference in total electronic energies, zero point vibrational energies and rotational and 
translational contributions at 298K. 

d) The CCSD values use the CISD zero point vibrational energies. 

e) Ref [70], 


Table 4 

Predicted binding energies (kcal/mol) and dipole moments (Debyes) 
for the anionic hydrogen bonded complexes. The binding energies 
were computed with respect to the most stable dissociation products 
as indicated in Table 3. 


Anionic 

complex 

Method 

Energy 

A E a 

AE b 

Dipole 

Moment 

FH CN- 

SCF 

-192.442703 

22.1 

22.0 

2.92 


MP2 

-193.056400 

27.8 

26.8 

- 


CISD 

-192.940474 

25.3 

24.3 

2.38 


CCSD' 

-193.059361 

25.5 

24.5 

- 

FH -NC- 

SCF 

-192.444070 

23.0 

22.9 

3.19 


MP2 

-193.054960 

26.1 

25.0 

- 


CISD 

-192.940721 

25.4 

24.3 

3.01 


CCSD' 

-193.059157 

25.4 

24.3 

- 

- 

Expt d 

- 

- 

21.1 

- 

H 2 0 CN- 

SCF 

-168.424382 

12.4 

12.3 

4.29 


MP2 

-169.036450 

16.3 

15.6 

- 


CISD 

-168.920387 

14.5 

13.8 

3.84 


CCSD' 

-169.043811 

15.2 

14.5 

- 

H 2 0 NC- 

SCF 

-168.425723 

13.4 

13.4 

4.48 


MP2 

-169.036487 

16.2 

15.4 

- 


CISD 

-168.921256 

15.1 

14.3 

4.32 


CCSD' 

-169.044430 

15.5 

14.7 

- 


Expt/ 



12.7±0.8 



a. Includes zero point energy and translational, rotational correction for 298K, see ref 71 for method. 

b. Includes zero point energy and translational, rotational correction and basis set superposition error 
determined by the counterpoise method, ref 53. 

c. CCSD energy performed at CCSD equilibrium geometry. The single point energy allowed all orbitals 
to be active, whereas in the geometry optimization the core and corresponding virtual orbitals were frozen. 
The optimum CCSD energies are *193.004605 and -193.004374 for FH CN” and FH NC”, respectively. 
CISD zero point energies were used. 

d. Ref 24. 

e. CCSD energy at the CISD equilibrium geometry. All orbitals active in the CCSD procedure. CISD zero 
point energies were used. 

f. Ref 29. 


Table 5 

Geometrical structures for the anionic hydrogen bonded complexes 
FH* "CN~ and FH ‘NC~. Bond lengths arc given in A. 


Anionic complex 

Method 

r Fir 

*CN 

TJT...C 

R-R..C 

FH-CN- 

SCF 

0.9472 

1.1457 

1.7790 

2.7262 


MP2 

1.0103 

1.1803 

1.6039 

2.6142 


CISD 

0.9769 

1.1585 

1.6768 

2.6537 


CCSD 

0.9879 

1.1696 

1.6666 

2.6546 

• 


*Fir 

re/v 

*n...N 

R F.../V 

a 

a 

o 

i 

SCF 

0.9437 

1.1483 

1.6481 

2.5918 


MP2 

0.9898 

1.1820 

1.5460 

2.5358 


CISD 

0.9678 

1.1611 

1.5798 

2.5476 


CCSD 0.9768 1.1719 1.5785 2.5553 



Table 6 

Geometrical structures for the anionic hydrogen bonded complexes 
H 2 0 CN“ and H 2 0 NC~. See Figure 2 for definitions of 
the molecular bond angles. Bond lengths are given in A 
and bond angles in degrees. 


Anionic 

complex 

Method 

r oh x 

r oh? 

rciv 

rc./f, 

Rco 

Or 

0 

7 

H 2 0 CN- 

SCF 

0.9599 

0.9392 

1.1484 

2.1246 

3.0706 

176.3° 

00 

to 

0 

103.6° 


MP2 

0.9966 

0.9586 

1.1836 

1.9055 

2.8984 

176.8° 

173.9° 

102.0° 


CISD 

0.9769 

0.9486 

1.1614 

1.9903 

2.9601 

176.6° 

171.5° 

O 

to 

-1 

0 



*OH x 

TOHa 

fCiV 

*N...H X 

R NO 

a 

0 

7 

H 2 0 NC" 

• SCF 

0.9596 

0.9390 

1.1500 

1.9501 

2.8963 

173.0° 

168.3° 

103.6° 


MP2 

0.9914 

0.9581 

1.1847 

1.8012 

2.7885 

h- » 

to 

bo 

0 

173.6° 

102.2° 


CISD 0.9749 0.9484 1.1629 1.8523 2.8199 172.3° 171.4° 102.8' 


Table 7 

Harmonic vibrational frequencies and infrared intensities 
for the anionic hydrogen bonded complexes. 

Frequencies are given in cm -1 and IR intensities (in parentheses) in km/molo . 


Anionic 

Complex 

Normal 

mode 

SCF 

MP 2 

CISD 

Expt. a 

FH-CN- 


3320 (2203) 

2352 

2844 ( 2679) 

1800 


) 

2373 

( 21 ) 

2044 

2261 

(7) 

2500 


Ws(<r) 

246 

(46) 

298 

278 

(61) 

- 


u 4 (it) 

1043 

(190) 

1142 

1107 

(143) 

1100 


■>(*■) 

162 

(3) 

159 

163 

(3) 

- 

FH-NC- 


3406 (2271) 

2727 

3034 (2618) 

1800 


CJ 2 ((7) 

2346 

( 101 ) 

2030 

2235 

(74) 

2500 



276 

(58) 

309 

300 

(70) 

- 


u 4 {ir) 

1034 

(254) 

1104 

1086 

( 211 ) 

1100 


w s (ir) 

123 

(5) 

127 

126 

(8) 

- 

h 2 o-cn- 


4195 

(43) 

3911 

4050 

(28) 



Ll? 2 ( ^ ) 

3786 

(801) 

3174 

3497 

(1106) 



Uji ( a/) 

2345 

(30) 

2012 

2232 

( 12 ) 



o,’ 4 (a/) 

1822 

(114) 

1722 

1777 

( 88 ) 



W.%(a/) 

406 

(73) 

472 

446 

( 66 ) 



o; fi (a/) 

175 

(26) 

217 

202 

(32) 



u >7 (a/) 

96 

(3) 

93 

98 

(3) 



u;fi(a//) 

798 

(116) 

902 

860 

(84) 



0 J 9 ( a/ / ) 

115 

(13) 

106 

114 

( 10 ) 


h 2 o-nc- 

u;i(a/) 

4198 

(43) 

3919 

4054 

(30) 



w 2 (a/) 

3800 

(864) 

3289 

3548 1 

(1116) 



UJ 3 ( a/ ) 

2329 

(74) 

2004 

2216 

(50) 



w 4 (a/) 

1827 

(127) 

1730 

1782 

( 100 ) 



w s (a/) 

398 

(73) 

463 

439 

( 66 ) 




197 

(34) 

237 

225 

(41) 



( a/) 

71 

(4) 

75 

72 

(5) 



w»(a//) 

802 

(133) 

890 

856 

( 101 ) 



u- f 9 (a// ) 

90 

(30) 

96 

93 

(33) 



a. Experimental fundamental frequencies are taken from ref [28]. 


Table 8 

Rotational Constants (MHz) for the Equilibrium structures of the 

Anionic Complexes. 


Anionic 

Complex 

SCF 

MP2 

CISD 

CCSD 

FH-CN" 

A 

3825 

4021 

3966 

3946 

FH-NC" 

A 

h 2 o-cn- 

4348 

4438 

4443 

4404 

' A 

513050 

536470 

530720 


B 

3370 

3631 

3548 


C 

H 2 0 - NC- 

3348 

3607 

3524 



A 

B 

C 


455610 

3887 

3854 


469440 

4065 

4031 


455490 

4033 

3998 



Figure Captions 


Figure 1. Definition of the internal coordinates for the H20'"CN anionic hydrogen- 
bonded complex. 

Figure 2. Definition of the internal coordinates for the H 2 0 "NC - anionic hydrogen- 
bonded complex. 

Figure 3. Valence electron density difference plot for the FH 'CN“ anionic complex. 
Short dashed lines indicate a depletion of electron density while solid fines indicate 
an increase of electron density. 

Figure 4. Valence electron density difference plot for the FH ' NC - anionic complex. 
Short dashed fines indicate a depletion of electron density while solid fines indicate 
an increase of electron density. 

Figure 5. Valence electron density difference plot for the HCN'HF hydrogen- 
bonded complex. Short dashed fines indicate a depletion of electron density while 
solid fines indicate an increase of electron density. 
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Abstract 

It was recently proposed that the Euclidian norm of the t 3 vector of the 
coupled cluster wave function (normalized by the number of electrons included in 
the correlation procedure) could be used to determine whether a single-reference- 
based electron correlation procedure is appropriate. This diagnostic, 7i, is defined 
for use with self-consistent-field molecular orbitals and is invariant to the same 
orbital rotations as the coupled cluster energy. T 3 is investigated for several different 
chemical systems which exhibit a range of multireference behavior, and is shown 
to be an excellent measure of the importance of non-dynamical electron correlation 
and is far superior to Co from a singles and doubles configuration interaction wave 
function. It is further suggested that when the aim is to recover a large fraction of 
the dynamical electron correlation energy, a large T\ (i.e., > 0.02) probably indicates 
the need for a multireference electron correlation procedure. 


1 Mailing Address: NASA Ames Research Center, Moffett Field, California 94035 
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Introduction 


It was recently proposed 1 that the Euclidian norm of the vector of tj ampli- 
tudes in the closed-shell coupled cluster singles and doubles wave function could be 
used as a diagnostic for the a priori prediction of the reliability of results obtained 
from a single-reference-based electron correlation procedure. The t 3 amplitudes in 
coupled cluster theory are closely related to the coefficients of singly excited config- 
urations in configuration interaction theory. It is well documented 2 that the singly 
excited configurations in an electron correlation procedure allow molecular orbital 
relaxation to occur. For many years quantum chemists have used Co, the reference 
configuration coefficient in a configuration interaction wave function, as a diagnos- 
tic. As is widely recognized, however, if Co is taken from a self-consistent-field 
(SCF) singles and doubles configuration interaction (CISD) wave function, then 
it is of limited utility since the molecular orbitals are strongly biased towards the 
SCF reference function. Thus it is not uncommon for a known multireference sys- 
tem to yield an SCF-CISD Co which is 0.95 or larger (i.e., the SCF determinant 
comprises 90% of the wave function). A reliable diagnostic which is more sensitive 
to the importance of non-dynamical electron correlation would therefore be of great 
utility. 

Laidig, Purvis and Bartlett 3 ' 4 have investigated the use of localized molecular 
orbitals in coupled cluster methods, specifically the doubles and the singles and 
doubles coupled cluster methods (CCD and CCSD, respectively). We note that 
the particular localization technique investigated by Laidig et al. does not leave 
the SCF energy unaffected 3,4 . They found that the use of localized molecular 
orbitals greatly improved the CCD results, but that the CCSD energies were little 
affected by the different reference molecular orbitals. The inclusion of e Tl in the 
CCSD wave function thus accounts for the important orbital relaxation effects which 
were incorporated by localizing the molecular orbitals. In addition, Scuseria and 
Schaefer 0 have investigated the use of Brueckner-like molecular orbitals in CCSD 
and CCD calculations and arrive at essentially the same conclusions. 

The purpose of the present study is to further investigate the use of the Euclid- 
ian norm of C as a diagnostic; applying this test to chemical systems exhibiting a 
range of bonding situations and known multireference and strongly single-reference 
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dominated problems. In this way the actual value and utility of the Euclidian norm 
of tj as a diagnostic tool will become evident. To begin, it must be emphasized 
that the diagnostic reported here 


was always determined using SCF molecular orbitals. As the results of references 
3 through 5 clearly demonstrate, it is possible to obtain a similar CCSD energy 
with different molecular orbitals which will give a different tj vector and a different 
Euclidian norm. In fact, the Euclidian norm of tj for the ‘optimized orbitals’ 
of reference 5 should be very close to zero. Thus, in order to compare T : from 
different chemical systems, the diagnostic must be uniquely defined for each system. 
The most straightforward .approach is to require that restricted Hartree-Fock SCF 
molecular orbitals are used to determine the CCSD wave function and diagnostic 
for each system, and this is therefore the approach which has been adopted in the 
present study. In addition, we point out that since the CCSD energy is invariant to 
unitary transformations of occupied-occupied or virtual-virtual molecular orbitals 
the Tj diagnostic will also be invariant to these types of orbital rotations. 

The next section contains a brief summary of the theoretical methods used 
together with a description of the method we have devised to judge the Tj diag- 
nostic. The results, including a discussion, are presented in the third section. Our 
conclusions are presented in the final section. 

Methods 

All of the chemical systems included in this study have been investigated 
previously 1,6 ' 7 and these reports include a detailed description of the basis sets 
and geometries. We therefore include only a brief description of the basis sets. 
Table 1 contains the size of the primitive basis, our designation and the reference 
from which the orbital exponents and contraction coefficients may be obtained. 
In forming the designation for each basis two rules have been followed. Firstly, a 
generally contracted atomic natural orbital (ANO) basis set is denoted by square 
brackets, e.g., [4321], where the numbers enumerate the number of contracted s, p, 
d and / functions, respectively. Secondly, a basis set which utilizes a segmented 
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contraction scheme is designated as 7s3p2dlf, for example. In most cases where a 
segmented contraction is used, the contraction has been performed over the core 
atomic orbitals, allowing maximum flexibility in the valence region. For those cases 
where the polarization function orbital exponents are not given in the reference 
the exponents are listed in Table 1. In addition, where more than one level of 
polarization function has been included (e.g., 7s3p2dlf Be) the levels are separated 
by a semicolon. 

Bond lengths are given in atomic units, a 0 . The unique bond length is specified 
for the Be 3 , Mg 3 , Be 4 and Mg 4 clusters. The trimers form an equilateral triangle 
and the tetramers adopt a tetrahedral structure. The pentamer, Bes, is defined by 
two bond lengths since it conforms to a trigonal bipyramidal geometry. For this 
system, the first bond length refers to a side of the triangular base while the second 
refers to the distance from an apex atom to one contained in the base. The bond 
lengths and bond angles for FOOF, (N0)2 and FNNF are the TZ2P MP2 structures 
reported in reference 1. 

Since the definition of T\ depends upon the number of electrons correlated 
it is clearly important to consider which electrons should be included in this def- 
inition. It is expected that only the valence electrons should be important for 
non-dynamical electron correlation effects and therefore we have chosen to freeze 
the core-like molecular orbitals in all procedures. It is possible that even if the core 
electrons are included in the correlation procedure, then the definition of T\ should 
include only the number of valence electrons (see note added in proof in reference 1). 
However, as several studies have demonstrated 2 * 8 ' 9 the basis set requirements for the 
adequate treatment of core- valence and core-core correlation effects are quite severe. 
Therefore, for our initial investigations of T\ only the valence electrons are consid- 
ered. Additionally, for basis sets which utilize segemented contractions the virtual 
molecular orbitals which are the core-counterparts were deleted from the correlation 
procedure. The CCSD wave functions were determined with a vectorized closed- 
shell CCSD method 10 and the Cl wave functions were evaluated with either the 
Berkeley shape-driven graphical unitary group Cl program 11 or the MOLECULE- 
SWEDEN codes 12 - 13 . 
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Results and Discussion 


The Ti diagnostic together with C 0 from CISD and full Cl wave functions are 
presented for several systems in Table 2. Note that only two electrons are correlated 
for the first five systems. Comparing T x with C 0 for these systems it is clear that 
there is a good correspondence betwreen and the total weight of the reference in 
the full Cl wave function. Thus for He and H 2 , where C 0 is greater than 0.99, the 
Tj diagnostic is 0.0029 and 0.0050, respectively, whereas for the other three systems 
(Be, Mg and Li 2 ) 7i is greater than 0.015 and C 0 is less than 0.965. Be and Mg are 
known to exhibit multireference behaviour due to the s — p near degeneracy. Li 2 
possesses a a — a* near degeneracy in addition to the s — p near degeneracy. 

Since the remaining molecules in Table 2 (He 2 , Be 2 , Mg 2 and HF) all contain 
more than two valence electrons it is possible to compare T 3 , the C 0 from CISD, and 
the C 0 from full Cl. As expected, He 2 is strongly dominated by a single reference 
function and so there is not a significant difference between the CISD and full Cl 
C 0 . Consistently, T 3 is again very small and is actually the same (to the precision 
reported) as for the single He atom. However, for the Be 2 , Mg 2 and HF diatomics 
there is a significant difference between the full Cl and CISD Co- In fact, for 
Be 2 the difference amounts to 4.5% of the full Cl wave function. An important 
point which should be emphasized is that due to the lack of size-extensivitv the 
discrepancy between the full Cl and CISD C 0 is expected to become larger as the 
number of electrons correlated increases. is greater than 0.013 for Be 2 , Mg 2 and 
HF demonstrating that a large degree of orbital relaxation occurs. Thus, the results 
of Table 2 demonstrate two important points: 1) there is a good correspondence 
between T 3 and the full Cl Co when a modest number of electrons are correlated 
and 2) for chemical systems with more than two electrons there may be a large 
difference between the CISD and full Cl Co. 

Tj and the CISD Co for several different chemical systems with a large number 
of valence electrons are collected in Table 3. The 7 " ] diagnostic and the Co for the 
Be and Mg clusters (at their equilibrium structures) indicates that these systems 
are probably not well described by a single-reference method and that a large degree 
of orbital relaxation is taking place. Binding energies and equilibrium bond lengths 
for the clusters, for example, would be expected to be substantially in error when 
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a single-reference-based treatment is used, and it is doubtful that binding energy 
predictions would be reliable to within even 10 kcal/mol. However, the large Co 
for these systems might tempt many observers to believe that a single-reference- 
based electron correlation procedure is adequate. Conversely, the Tj value is larger 
than 0.02 for each cluster with the exception of Mg 3 . The comparisons made in 
Table 2 together with multireference Cl (MRCI) results 6 suggest that multireference 
techniques are required for Be and Mg clusters, and thus that a value larger 
than 0.02 is a clear indication that other important configurations exist and may be 
needed as references m a treatment of dynamical electron correlation. The infinite 
separation results for the Be and Mg clusters also demonstrate the inadequacy of 
using Co from a CISD w r ave function since the Co suggests that as the number 
of atoms increases the 'super-molecule 5 is more difficult to describe whereas the 
size-extensive CCSD method correctly show's that these systems are all equivalently 
described (in fact, since only valence electrons are correlated the CCSD results 
correspond to a full Cl). 

The FOOF, (NO) 2 and FNNF molecules are included in Table 3 since these 
were the systems investigated in the study 1 w'hich first suggested the use of Tj . 
These systems are very difficult to describe — the geometry of FOOF is not even 
qualitatively correct at the CISD level, for example. The C 0 values for FOOF and 
(NO), are very similar although Tj indicates that non-dynamical electron corre- 
lation is much more important for FOOF. The results of several single-reference 
methods for these twm systems 1 provide additional evidence that the electron cor- 
relation of FOOF is indeed even more difficult to describe than that for (NO) 2 . 
The Co for the isomers of FNNF suggests that these systems are more strongly 
dominated by a single reference than either FOOF or (NO) 2 and that they are all 
nearly equally well described by a single-reference-based method. However, wffiile 
the Tj diagnostic does suggest that non-dynamical electron correlation is less im- 
portant in the cis and trails isomers, it also indicates that the transition state is 
strongly affected by non-dynamical electron correlation and thus, single-reference- 
based methods will not w 7 ork as well for TS-FNNF as they do for cis and trans 
FNNF. Again, the latter conclusions are consistent with the results of reference 1. 

The last four molecules of Table 3 are all known to be strongly dominated 
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by a single determinant reference function and both Co and Tj are consistent with 
this observation. However, the fine details of relating Co and T] exhibit small 
inconsistencies. For example, for the first-row closed-shell hydrides it is generally 
accepted that the reliability of a single-reference-based electron correlation method 
decreases in the order CH 4 > H 2 0 > HF. The T a diagnostic is consistent with this 
empirical observation whereas the C 0 from a CISD wave function exhibits exactly 
the opposite trend. 

Perhaps the molecule which best exhibits the superiority of 7i over C 0 as a 
diagnostic is the CuH diatomic. It has been shown 14 that the bonding of CuH is 
complicated because of the importance of both the and (f 10 ^ 1 atomic occu- 

pations of Cu. Thus there are several important configurations which differ from 
the closed shell single determinant reference by a single .excitation. The C 0 for 
this diatomic is 0.96, which is very similar to that obtained for CH 4 . However, 
the Tj value, 0.046, is the largest found in this study. Thus, the important non- 
dynamical electron correlation effects present in the bonding of CuH are completely 
missed by the single-reference CISD method whereas the Tj diagnostic correctly 
indicates the importance of these effects. 

Table 4 contains Tj and C 0 (from CISD) for Be 3 and H 2 0 using several basis 
sets in order to determine the one-particle basis set effect. In order for this diagnostic 
to be generally useful it should exhibit a certain degree of invariance with respect 
to the choice of a one-particle basis set. This statement assumes, of course, that 
the smallest basis sets at least contain proper correlating functions. On the other 
hand, it is well known that the one-particle and n-particle basis sets are inherently 
coupled though this coupling is usually small. As the one-particle basis set limit 
is approached, it may be expected that T 2 will stabilize. This is expected despite 
the fact that the n-particle basis increases substantially with one-particle basis set 
augmentations. 

The results of Table 4 confirm the above discussion and demonstrate that 
Ti converges to a value near 0.0340 for Be 3 and near 0.0075 for H 2 0. The fact 
that Tj decreases with improvements in the one-particle basis set provides further 
support for the above discussion. In other words, a larger degree of orbit al relaxation 
is required for the smaller one-particle basis sets (giving a larger Tj ) in order to 


compensate for the lack of flexibility. Thus the value to which X converges should 
give an indication of the inherent importance of non-dvnamical electron correlation 
for the chemical system under investigation. Moreover, the rate of convergence with 
respect to basis set improvement should give a measure of the interaction between 
the one- and n-particle basis sets. 


Conclusions 

The X] diagnostic has been shown to be a reliable measure of the impor- 
tance of non-dynamical electron correlation and to be far superior to the use of C 0 
from a CISD wave function as an indicator as to whether it is appropriate to use 
a single-reference-based electron correlation procedure. Iko doubt a similar type of 
diagnostic could be defined for the CISD wave function by separating the Cj coeffi- 
cients (coefficients from the singly excited configurations). However, there are two 
problems with this procedure. Because CISD is, in general, not size extensive then 
the diagnostic would not have the desired property of giving the same result for two 
non-interacting He atoms as it would for a single He atom. Also, the results for 
CuH presented in this study clearly indicate that the CISD procedure is incapable 
of overcoming the bias of using SCF molecular orbitals and thus any diagnostic sim- 
ilar to X] but. based on an SCF-CISD wave function would almost certainly suffer 
from this bias. 

Several studies 14,22 have pointed out that the coupled pair functional 23 (CPF), 
modified CPF 24 (MCPF) and averaged CPF 25 (ACPF) methods are very good at 
identifying specific configurations which are important and hence should be used as 
references in a multireference electron correlation procedure. A diagnostic similar 
to X i could also be constructed for these methods and it is likely that it would give 
similar results to X] for those situations where the CCbD and the various CPF- 
tvpe methods gave similar results. Clearly for the situation where the methods 
give very different results (such as FOOF 1 ) the diagnostics would be expected 
to yield different results also. In addition these observations also indicate that 
specific important configurations may be identified by analysis of the CCSD D and 
X amplitudes. 

Finally, the results of this study indicate that if Xj is greater than 0.02 then 
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single-reference-based electron correlation methods are probably unreliable and will 
certainly not yield highly accurate results. 
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Table 1 

Basis set designations and definitions used in this study. 


Atom 

Primitive Basis 

Designation 

Reference 

Polarization 

Exponents 

H 

4s 

DZP 

15,16 

0.75 

H 

8s2p 

6s2p 

17 

1.0,0.33 

H 

8s6p 

[32] 

18 

- 

H 

8s6p4d 

[321] 

18 

- 

H 

8s6p4d 

. [432] 

18 

- 

He 

8s2p 

6s2p 

17 

1.0,0.33 

Li 

9s4p 

4s3p 

17 

- 

Be 

12s5p2d 

7s3p2d 

17 

0.3, 0.1 

Be 

12s5p2dlf 

7s3p2dlf 

17 

0.3,0.1;0.26 

Be 

12s7p4d2f 

[5321] 

18 

- 

Be 

12s7p4d2f 

[6432] 

18 

- 

C 

10s6p2d 

TZ2P Q 

15,19 

1.5,0.35 

N 

10s6p2d 

TZ2P“ 

15,19 

1 .5,0.35 

0 

10s6p2d 

TZ2P° 

15,19 

1.5,0.35 

0 

13s8p6d 

[432] 

18 

- 

0 

13s8p6d4f 

[4321] 

18 

- 

0 

13s8p6d4f2g 

[54321] 

18 

- 

F 

9s5pld 

DZP 

15,16 

1.6 

F 

10s6p2d 

TZ2P a 

15,19 

1.5.0.35 

Ne 

10s6p2d 

TZ2P 6 

17 

4.5, 1.3 

Mg 

12s9p2d 

6s5p2d 

20 

0.3, 0.1 

Cu 

14sllp6d3f 

8s6p4dlf 

7 



a. The 5s3p contraction of reference 19 was used. 

b. A 5s3p contraction, similar to those given in reference 19, was constructed. 
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Table 2 

The Tj diagnostic together with the Co obtained 
from a full Cl and a CISD wave function 0 . 


Molecule 

Basis Set 

T 

T, 

c b 

cs 

He 

6s2p 

- 

0.0029 

0.9960 

0.9960 

h 2 

6s2p 

1.361 

0.0050 

0.9912 

0.9912 

Be 

7s3p2d 

- 

0.0210 

0.9523 

0.9523 

Mg 

6s5p2d 

- 

0.0159 

0.9640 

0.9640 

Li 2 

4s3p 

5.11 

0.0165 

0.9510 

0.9510 

He 2 

6s2p 

5.61 

0.0029 

0.9920 

0.9921 

Be 2 

7s3p2d 

4.75 

0.0282 

0.8901 

0.9150 

Mg 2 

6s5p2d 

7.35 

0.0138 

0.9268 

0.9401 

HF rf 

DZP 

2.5995 

0.0187 

0.9583 

0.9680 


a. All correlated wave functions are based upon SCF molecular orbitals. Only 
valence electrons have been included in the correlation procedure. Bond lengths are 
in atomic units. Go. 

b. FCI. 

c. CISD. 

d. FCI and CISD results from reference 21. 
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Table 3 


The Tj diagnostic together with C 0 obtained from a CISD wave function 3 . 


Molecule 

Basis Set 

Geometry 

T, 

Co 

CuH 

b 

2.850 

0.0461 

0.9621 

Be 3 

7s3p2d 

4.273 

0.0360 

0.9133 

Be 4 

7s3p2d 

3.915 

0.0318 

0.9189 

Be 5 

7s3p2d 

3.831,3.929 

0.0290 

0.9094 

Be 3 

7s3p2d 

OG 

0.0210 

0.9067 

Be 4 

7s3p2d 

CO 

0.0210 

0.8933 

Be 3 

7s3p2d 

OO 

0.0210 

0.8828 

Mgs 

6s5p2d 

7.522 

0.0127 

0.9235 

Mg4 

6s5p2d 

6.102 

0.0204 

0.9102 

Mg 3 

6s5p2d 

OO 

0.0159 

0.9240 

Mg4 

6s5p2d 

OO 

0.0159 

0.9111 

FOOF 

TZ2P 

see text 

0.0313 

0.9189 

(NO) 2 

TZ2P 

see text. 

0.0203 

0.9177 

cis-FNNF 

TZ2P 

see text 

0.0187 

0.9303 

trans-FNKF 

TZ2P 

see text 

0.0166 

0.9308 

TS-FNNF C 

TZ2P 

see t ext. 

0.0277 

0.9283 

HF 

TZ2P 

1.734 

0.0104 

0.9775 

H 2 0 

TZ2P 

1.809.104.8° 

0.0096 

0.9720 

ch 4 

TZ2P 

2.052 

0.0073 

0.9672 

Ne 

TZ2P 

- 

0.0065 

0.9850 


a. All correlated wave functions are based upon SCF molecular orbitals. Onlv 
valence electrons have been included in the correlation procedure. Bond lengths are 
in atomic units, a 0 . 

b. The Cu basis is as described in Table 1 and the H basis is the [32] ANO basis 
set. 

c. Transition state to cis-trans isomerization. 
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Table 4 

The 7i diagnostic for Be 3 and H 2 O using several different 
basis sets. The Cq value is obtained from a CISD wave function. 


Molecule 

Basis Set 

T, 

Co 

Be 3 

7s3p2d 

0.0360 

0.9133 

Be 3 

7s3p2dlf 

0.0339 

0.9149 

Be 3 

[421] 

0.0386 

0.9107 

Be 3 

[5321] 

0.0341 

0.9148 

Be 3 

[6432] 

0.0341 

0.9157 

H 2 0 

TZ2P 

0.0096 

0.9720 

h 2 o 

[432/32] 

0.0076 

0.9721 

h 2 o 

[4321/321] 

0.0071 

0.9714 

h 2 o 

[54321/432] 

0.0078 

0.9713 


a. The geometries are the same as those listed in Table 2. Only valence electrons 
have been correlated. 
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